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Abstract:  In  this  paper,  a  hysteresis  model-based 
nonlinear  robust  controller  is  developed  for  a  piezo¬ 
electric  actuator,  utilizing  a  Lyapunov-based  stability 
analysis,  which  ensures  that  a  desired  displacement  tra¬ 
jectory  is  accurately  tracked.  Simulation  results  are 
presented  and  discussed  to  demonstrate  the  proof  of 
concept  for  the  proposed  robust  control  strategy. 

1  Introduction 

Piezoelectric  actuator  (PZTA)  based  systems  are 
emerging  as  an  important  technology  for  precise  po¬ 
sitioning  and  have  received  wide  attention  in  both 
the  scientific  and  industrial  communities.  These  de¬ 
vices  are  capable  of  completing  high  precision  actua¬ 
tion  tasks.  They  are  often  utilized  in  motion  actuation 
applications,  due  to  their  high  stiffness,  fast  response, 
and  physically  unlimited  resolution  [19],  and  can  be 
used  as  either  sensors  or  actuators  in  control  systems. 
The  advantages  of  the  PZTA  [14]  include:  1)  no  wear, 
2)  high  efficiency,  3)  almost  infinite  small  positioning 
ability,  4)  ultra  fast  expansion,  and  5)  capability  to 
deliver  large  actuation  forces.  These  micro-positioning 
elements  play  a  big  role  in  many  vital  applications  such 
as  scanning  tunneling  microscopy  [31],  scanning  probe 
microscopy  ([9]  and  [23]),  laser  applications  [29],  and 
hydraulic  servo  control  systems  [5]. 

Despite  the  advantages  of  PZTA’s  for  micro  position¬ 
ing  applications,  the  actuator’s  positioning  response 
shows  a  strong  hysterical  activity  due  to  their  composi¬ 
tion  from  ferroelectric  ceramic  materials.  Specifically, 
an  applied  voltage  is  typically  the  input  control  signal 
which  activates  the  PZTA.  In  the  event  that  the  input 
control  voltage  is  relatively  large,  the  PZTA  exhibits 
a  significant  amount  of  distortion  due  to  the  inherent 
hysteresis  in  the  device,  and  this  effect  may  reduce  the 
stability  of  the  system  in  feedback  control  applications 
[26].  Due  to  this  nonlinear  behavior,  one  would  expect 
difficulties  in  using  PZTA’s  for  precise  tracking  control 
applications.  Hence,  nonlinear  control  strategies  are 
needed  for  the  use  of  the  PZTA’s  in  micro  positioning 
and  tracking  systems. 
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Past  PZTA  research  has  focused  on  establishing  an  ac¬ 
curate  dynamic  actuator  model  ([3],  [18],  and  [30]), 
while  other  research  has  focused  on  the  development 
of  active  control  strategies  for  use  in  precision  posi¬ 
tioning  and  tracking  applications  ([13],  [14],  [19],  [25], 
and  [27]).  In  [12],  the  authors  provided  a  concise  lit¬ 
erature  summery  with  regard  to  modeling  the  PZTA’s 
dynamics.  Recently,  the  literature  has  focused  on  the 
development  of  alternate  models  to  describe  the  hys¬ 
teresis  within  the  PZTA,  due  to  the  changeling  non¬ 
linear  nature  of  this  phenomenon.  In  [30],  the  authors 
included  a  nonlinear  spring  element  into  the  hysteresis 
model  and  utilized  a  Maxwell-slip  structure,  while  the 
authors  of  [16]  used  a  support  vector  regression  non¬ 
linear  model  and  neural  networks.  The  authors  of  [17] 
and  [22]  found  that  voltage  to  displacement  lineariza¬ 
tion  of  a  PZTA  may  be  achieved  if  the  control  input  is 
the  applied  electric  charge  rather  than  the  applied  volt¬ 
age.  Furutani  et  al.  [10]  was  able  to  improve  the  PZTA 
control  strategy  by  combining  induced  charge  feedback 
with  inverse  transfer  function  compensation.  Vautier 
and  Moheimani  [26]  showed  that  applying  an  electric 
charge  to  control  the  position  reduces  the  effects  of  the 
nonlinearity  in  the  PZTA  dynamics.  They  also  showed 
experimental  results  demonstrating  the  effectiveness  of 
using  electrical  charge.  Further,  Main  et.  al.  [21]  pre¬ 
sented  experimental  data  from  tests  for  both  voltage 
and  charge  control  which  showed  that  charge  control  is 
significantly  more  linear  and  less  hysteretic  than  volt¬ 
age  control  over  the  same  actuator  displacement  range. 

In  [20],  the  authors  proposed  a  new  mathematical 
model  to  describe  complex  hysteresis  that  is  based  on 
a  new  parameter  called  turning  voltage  of  a  PZTA.  In 
their  work,  the  authors  of  [20]  were  able  to  utilize  this 
parameter  to  suppress  the  inherent  hysteresis  to  within 
±1%  full  span  range  of  a  PZTA.  Shieh  et  al.  in  [24] 
developed  a  parametrized  hysteretic  friction  function, 
based  on  the  LuGre  model,  to  describe  the  PZTA  hys¬ 
teresis  behavior.  Thus,  these  researchers  were  able  to 
design  an  adaptive  displacement  tracking  control  with 
a  parameter  adaptation  algorithm.  In  [3],  Bashash  and 
Jalili  presented  a  perturbation  estimation  technique  to 
compensate  for  the  structural  nonlinearities  and  un¬ 
modeled  PZTA  dynamics.  In  their  work,  the  authors 
of  [3]  validated  the  proposed  model  by  experimental 
tests  using  a  PZTA-driven  nano-stager  with  capacitive 


position  sensor. 

Other  research  in  the  literature  has  focused  on  the  de¬ 
velopment  of  intelligent  control  schemes  for  the  preci¬ 
sion  control  of  PZTA’s.  Some  of  these  designs  were 
based  on  the  inverse  hysteresis  model  that  is  assumed 
to  be  known  a  priori,  so  feedforward  techniques  can 
be  utilized  in  the  control  design  ([2]  and  [18]).  Others, 
such  as  [4],  [6],  [13]  and  [14],  applied  a  feedback  lin¬ 
earization  to  compensate  for  the  hysteresis  dynamics 
and  then  a  tracking  controller  was  implemented.  Wu 
and  Zou  in  [28]  presented  an  inversion-based  iterative 
control  approach  to  compensate  for  both  the  hysteresis 
and  the  vibrational  dynamics  variations  during  high¬ 
speed,  large  range  tracking.  Neural  networks  and  fuzzy 
controls  were  also  utilized  to  model  the  PZTA  hystere¬ 
sis  nonlinearities  and  control  the  micro  motion  of  the 
PZTA  ([14]  and  [27]).  In  [19],  the  authors  proposed 
a  robust  control  strategy  for  precise  positioning  track¬ 
ing.  The  implementation  of  their  control  law  requires 
only  a  knowledge  of  the  estimated  system  parameters 
and  their  corresponding  bounds  as  well  as  the  bound 
of  the  hysteresis  effect  including  disturbances.  In  [25], 
Stepanenko  et  al.  introduced  and  implemented  an  ap¬ 
proximation  function  to  compensate  for  the  hysteresis 
nonlinearities  by  fuzzy  logic  techniques. 

In  this  paper,  the  displacement  of  a  PZTA  is  actively 
controlled  to  track  a  desired  trajectory.  A  nonlinear 
robust  control  strategy  is  developed  based  on  the  feed¬ 
back  of  the  PZTA  charge,  and  partial  knowledge  of  the 
hysteresis  model.  Charge  steering  approach  eliminates 
the  effect  of  the  hysteresis  in  the  PZTA  and  provides 
better  robust  control  over  the  voltage  steering  approach 
[1].  The  charge  measurement  is  obtained  by  measuring 
the  voltage  across  a  capacitor  that  is  added  in-series  to 
the  PZTA  circuit.  A  Lyapunov-based  analysis,  which 
proves  precise  tracking,  is  utilized  to  develop  the  con¬ 
trol  strategy.  Simulation  results  are  presented  which 
demonstrates  the  proof  of  concept  of  the  active  con¬ 
trol  approach.  This  paper  is  organized  as  follows.  In 
Section  2,  a  Coleman-Hodgdon-based  hysteresis  model 
along  with  the  piezoelectric  actuator  dynamics  are  pre¬ 
sented  as  well  as  the  required  assumptions  for  the  sys¬ 
tem.  In  Section  3,  a  nonlinear  robust  control  scheme  is 
developed  along  with  the  stability  analysis  which  veri¬ 
fies  that  the  piezoelectric  desired  displacement  can  be 
tracked.  In  Section  4,  simulation  results  are  presented 
for  the  robust  control  strategy.  Finally,  concluding  re¬ 
marks  are  provided  in  Section  5. 

2  PZTA  System  Model 

2.1  PZTA  Elongation  Dynamics 

A  PZTA  with  a  single  elongation  axis,  depicted  in  Fig¬ 
ure  1,  can  be  dynamically  modeled  as 

my  +  FL  =  Fp 


where  m  £  R  denotes  the  PZTA  mass,  L  £  R  de¬ 
notes  the  non-activated  length  of  the  PZTA,  Fp  ( t )  £  R 
denotes  the  force  generated  by  the  PZTA  elongation, 
Fl  (y,  y)  £  R  denotes  the  perpendicular  forces  acting 
on  the  PZTA,  y(t),y  (f)  ,y(t)  £  R  are  the  displace¬ 
ment,  velocity,  and  acceleration,  respectively,  of  the 
PZTA  effective  tip  of  elongation.  An  equivalent  cir¬ 
cuit  model  of  the  PZTA  can  be  described  as  shown  in 
Figure  2. 


Figure  1:  PZTA  elongation  induced  by  applied  voltage 


Figure  2:  PZTA  equivalent  circuit  model 

In  this  schematic,  V  (t)  £  R  denotes  the  applied  in¬ 
put  control  voltage,  q(t)  £  R  represents  the  current 
flowing  through  the  PZTA,  since  i(t)  =  q  (t) .  Further, 
Cm  £  R+  is  a  series  connected  capacitor  that  facili¬ 
tates  the  measurement  of  q  (t)  £  R,  through  a  mea¬ 
surement  of  the  voltage  Vm  ( t )  £  R  across  Cm,  recall 

p.  2 


(1) 


that  q  =  CmVm.  In  Figure  2,  the  circuit  element  indi¬ 
cated  by  H  models  the  inherent  hysteresis  between  the 
voltage  Vh  ( t )  £  R  and  the  induced  charge  q  (t).  A  sub¬ 
sequent  section  will  further  define  the  hysteresis  model. 
The  parameter  Cc  £  R+  is  the  internal  capacitance  of 
the  PZTA,  qc  ( t )  £  R  is  the  current  flowing  through 
Cc,  and  Vc  (t)  £  R  is  the  voltage  across  this  capaci¬ 
tance.  Finally,  the  circuit  element  indicated  by  Tem 
represents  the  subsequently  defined  elongation  model, 
where  Tem  £  R+  is  the  elongation  constant  [19],  and 
qp  ( t )  €  R  denotes  the  current  flowing  through  this  cir¬ 
cuit  branch. 

2.2  Hysteresis  Model 

A  nonlinear  hysteresis  model  H  (Vh)  £  R  can  be  defined 
to  describe  the  relationship  between  the  input  voltage 
Vh  ( t )  and  the  induced  charge  q(t).  For  this  paper,  a 
Duhem-based  hysteresis  model  can  be  defined  as  [8] 

q  =  H(Vh)±f(Vh)  +  d(Vh)  (2) 

where  /  (14)  £  R  is  a  subsequently  defined  signal.  The 
variable  d(Vh)  £  R  is  defined  as 

dine  =  [q0~  f{Vho)\e~5(y^o)  (3) 

vh 

+e  dVh  J  ( 9  (T)  _  (r))  e&TdT 

vho 

if  (vh  (i t )  >  0  and  Vh  (t)  >  Vho"j 

ddec  =  [qo~  f(Vho)}e-5(Vh°-Vh)  (4) 

vh 

+e5Vh  J  ( g  (t)  -  /  (r))  e~STdr 

vho 

if  (Vh  ( t )  <  0  and  Vh  (t)  <  Vh „) 

where  q0  £  R  is  the  induced  charge  at  t  =  t0,  Vha  £  R 
is  the  input  voltage  at  t  =  t0,  S  £  R+  is  a  con¬ 
stant,  and  /  (t)  £  R  is  the  partial  derivative  of  /  (14) 
with  respect  to  the  input  voltage  14  (t).  The  signals 
d  (14) ,  /  (14) ,  and  g  (14)  £  R  have  the  following  prop¬ 
erties  [7]: 

Property  1:  The  function  /  (14)  is  piecewise  smooth, 
monotone  increasing,  and  odd. 

Property  2:  The  function  g  (14  )  is  piecewise  contin¬ 
uous,  and  even. 

Property  3:  The  function  /  (14)  is  known,  and  in¬ 
vertible,  such  that  Vh  =  f  (Vh),  and  14  = 

r1 04)  • 

Property  4:  The  partial  derivative  of  /  (14)  with  re¬ 
spect  to  14(f)  ( i.e .,  d*Qy^)  is  not  identically  zero, 
hence,  /  (oo)  =  g0,  where  r/0  £  R+  is  a  constant. 


Property  5:  The  function  g  (Vh)  has  a  finite  up¬ 
per  limit  ^ i.e.,f  (oo)  =  g(oo)^  where  /  (14)  > 

9  (Vh). 

Property  6:  The  function  d  (Vh)  is  bounded  as 

Ml  <  Vd  (5) 

where  r]d  £  R+  is  a  constant  (see  Appendix  1  for 
proof). 

2.3  Elongation  Model 

Two  linear  relationships  model  the  elongation  action 
of  the  PZTA.  One  relationship  models  the  effect  of  the 
induced  charge  qp  (t)  on  the  displacement  of  the  elon¬ 
gation  axis,  y  (t) ,  and  can  be  written  as  [19] 

qP  =  Temy.  (6) 

The  second  relationship  depicts  the  force  Fp  (t)  im¬ 
parted  by  the  elongation  action  as  a  function  of  the 
voltage  Vc  (t)  which  becomes  [19] 

Fp  4  Tem Vc.  (7) 

In  (6)  and  (7),  Tem  is  the  elongation  constant  inherent 
to  the  PZTA. 

2.4  Dynamic  Model  and  Assumptions 

To  facilitate  the  subsequent  control  objective,  a  dy¬ 
namic  expression  is  desired  that  relates  the  displace¬ 
ment  of  the  elongation  axis  y  (t)  as  a  function  of  charge 
q  (t)  induced  within  the  PZTA.  The  advantage  of  work¬ 
ing  with  such  a  dynamic  model  is  clear  due  to  the  lack 
of  the  hysteresis  terms  which  has  been  discussed  in  the 
literature  (see  [10],  [21],  and  [26]).  From  Figure  2,  it  is 
clear  that  the  induced  charge  q  (t)  can  be  described  as 

q  =  CcVc  +  qp.  (8) 

From  (8),  the  expression  in  (1)  can  be  written  as 

my  +  FL=(J^jq  (9) 

where  (6)  and  (7)  were  utilized.  The  variable 
Fl  (y,y)  £  R  is  defined  as 

~Fl  —  FL  +  y.  (10) 

To  facilitate  the  tracking  control  design,  three  assump¬ 
tions  frame  the  analysis. 

Assumption  1:  The  PZTA’s  parameters  and  other 
parameters  to,  Cc,  Cm  and  Tem  are  assumed  to 
be  known,  and  constants  with  respect  to  time. 


Assumption  2:  The  velocity  y  (£)  and  displacement 
y(t)  of  the  PZTA  effective  tip  are  assumed  to  be 
measurable. 

Assumption  3:  It  is  assumed  that  the  forces  Fl  (y,  y) 
and  their  first  time  derivative  FL  (y,  y,  y)  are 
bounded  (i.e.,  FL  (y,  y) ,  FL  (y,  y,  y)  G  pro¬ 
vided  y(t),  y(t),y(t)  G  Too)- 

3  Robust  Control  Development 

3.1  Control  Design  Objective 

The  control  objective  is  to  ensure  that  the  displace¬ 
ment  y(t)  of  the  PZTA  effective  tip  tracks  the  desired 
trajectory  yd(t)  G  R  in  the  following  sense 

\yd(t)  -  y{t)\  <  £  as  t  y  oo  (11) 

where  e  G  R+  is  a  constant  that  can  be  selected  arbi¬ 
trary  small. 

Assumption  4:  The  subsequent  analysis  requires 
that  a  desired  trajectory  is  selected  such  that 
Vd(t),  Vd(t),  and  yd(t)  G  Loo  where  yd(t),yd(t) 
G  R. 

To  facilitate  the  subsequent  development,  a  filtered 
tracking  error  signal,  denoted  by  r(£)  G  R,  is  defined  as 

r  =  e  +  ae  (12) 

where  a  G  R+  is  a  control  gain.  The  variable  eft)  G  R 
is  defined  as 

e  =  yd-y-  (13) 

Based  on  the  definition  of  eft)  in  (13),  it  is  clear  that 
if  |e(£)|  <  £  as  £  — >•  oo,  then  | yd{t)  —  y(t)\  <  e  as 
£  — >  oo,  thus,  meeting  the  control  objective. 

3.2  Closed-Loop  Error  System 

To  facilitate  the  development  of  a  closed-loop  error  sys¬ 
tem,  a  control  strategy  must  be  developed  to  account 
for  the  inherent  hysteresis  that  exists  between  the  input 
voltage  Vh  (£)  and  the  induced  charge  q  (£).  To  continue 
this  development,  the  charge  as  defined  in  (2)  can  be 
rewritten  as 

q  =  f(V-Vm-Vc)  +  d(Vh).  (14) 

It  is  clear  from  Figure  2  that  an  expression  for  Vh  (£) 
can  be  stated  as 

vh  =  V  -  Vm  -  Vc.  (15) 

To  meet  the  previously  stated  control  objective,  the 
control  input  V  (£)  as  shown  in  Figure  2  can  be  designed 
as 

V  AVm+Vc  +  r1  (Vh)  (16) 

where  Vh  (£)  G  R  is  a  subsequently  designed  auxiliary 
control  signal,  and  Property  3  was  utilized. 


Remark  1  From  (8),  the  voltage  Vc  (£)  can  be  obtained 
as 

Vc  =  ^~{q-qP)  (17) 

where  the  charge  qp(t)  is  computed  from  (6),  and 
the  charge  q  (£)  is  computed  from  the  measurement  of 
Vm  (£)  across  the  capacitor  Cm  (i.e.,  q  =  CmVm). 

Utilizing  (16),  the  hysteresis  between  the  voltage  and 
change  in  (14)  can  be  written  as 

q  =  H(Vh)  =  Vh+d(Vh).  (18) 

To  continue  the  closed- loop  error  system  development, 
an  auxiliary  error  signal  s  (£)  G  R  is  defined  as 

s  =  qd~q  (19) 

where  qd  (y,  y)  G  R  is  a  subsequently  designed  auxil¬ 

iary  signal.  From  (18),  the  auxiliary  signal  s  (£)  can  be 
rewritten  as 

S  =  qd-Vh-  d(\ 4).  (20) 

The  auxiliary  control  signal,  V h  (£)  is  designed  as 

Vh  =  qd  +  kss  (21) 

where  ks  G  R+  is  a  constant. 

Remark  2  The  expression  in  (21)  can  be  expanded  as 

Vh  =  qd  +  ks(qd-q)  (22) 

where  (19)  is  utilized,  and  the  charge  q(t)  is  computed 
from  the  measurement  of  Vm  (£)  across  the  capacitor 
Cm  (i.e.,  q  =  CmVm) • 

From  (21),  the  expression  in  (20)  can  be  simplified  as 

s  =  -kss  -  d(Vh) .  (23) 

The  auxiliary  error  signal  s  (£)  can  be  upper  bounded 
as 

(24) 

where  r/s  G  R+  is  a  constant  and  Property  6  was  uti¬ 
lized. 

With  the  control  concept  in  place  to  account  for  the 
voltage  to  charge  hysteresis,  the  dynamic  model  as  de¬ 
fined  in  (9)  is  now  incorporated  to  complete  the  closed- 
loop  error  system  development.  From  the  PZTA  dy¬ 
namics,  described  in  (9),  the  following  expression  can 
be  obtained 

“S+Mlc)'  =  (if)'w  (25) 


where  (19)  was  utilized.  To  facilitate  the  development 
of  the  closed-loop  error  system  for  r(t),  the  first  time 
derivative  of  (12)  is  taken  and  then  both  sides  are  mul¬ 
tiplied  by  to,  thus,  obtaining  the  following  expression 

mr  =  myd  +  FL+  s  -  2a  +  am &  (26) 

where  the  first  time  derivative  of  (13),  and  (25)  were 
utilized.  To  facilitate  the  subsequent  analysis,  the  ex¬ 
pression  in  (26)  is  rewritten  as 

mr  =N  +Nd-e+  (^j  s  (^j  qd  (27) 

where  the  auxiliary  signal  N  ( y ,  y)  €  R  is  defined  as 

N=N-  Nd.  (28) 

The  variable  N(y,  y,  yd )  e  R  is  defined  as 

N  =  myd  +  F  l  +  e  +  ame  (29) 

and  Nd(t)  €  R  is  defined  as 

Nd  =  N\y=ydiy=yd  (30) 

=  myd  +  FLd 

where  F Ld  (yd,  yd )  €  R,  introduced  in  (10),  is  evaluated 
at  yd  and  yd.  Based  on  (27),  the  signal  qd(y,y )  is 
designed  as 


where  Ay  €  R+  is  a  constant  gain,  e  £  R+  is  a  small 
constant,  and  p(||z||)  €  R  is  a  function  of  norm  z  (t)  € 
R2.  The  variable  2  (t)  is  defined  as 

z  =  [  e  r  ]T  (32) 

where  r  (t)  and  eft)  were  introduced  in  (12)  and  (13), 
respectively.  By  utilizing  (28)- (30),  the  following  in¬ 
equality  can  be  developed  (see  Appendix  3  for  further 
details) 

^j<p(NI)NI-  (33) 

After  substituting  (31)  into  (27),  the  following  closed- 
loop  error  system  can  be  obtained 

mr  =  N  +Nd  -  e  +  s  -  krr  (34) 

-\p{  INI)2  INI2  r- 

3.3  Stability  Analysis 

Theorem  1  The  controller  given  in  (16),  (21),  and 
(31)  ensures  that  \e(t)\  <  e  as  t  — >  oo  provided  the 
control  gain  Ay,  introduced  in  (31),  is  sufficiently  large, 
hence,  e  ( t )  is  practically  regulated  to  zero.  A  secondary 
objective  is  that  all  closed-loop  signals  are  also  bounded 

Proof:  See  Appendix  2. 


4  Simulation  Results 

A  numerical  simulation  was  performed  to  demonstrate 
the  proof  of  concept  of  the  proposed  robust  tracking 
controller.  The  PZTA  dynamic  model,  introduced  in 
(25),  was  used  to  implement  the  actuator  in  the  sim¬ 
ulation.  The  simulated  PZTA  parameters  were  set 
equal  to  the  following:  to  =  0.01  [kg],  Cc  =  1.5  [yF], 
and  Tem  =  1  [/, iN/V].  The  measurement  capacitor  Cm 
was  omitted  for  this  simulation  due  to  ability  to  mea¬ 
sure  the  charge  q  (t)  within  the  simulation  environment. 
The  forces  FL  (y,  y)  on  the  PZTA  was  defined  as 

Fl  =  ky  +  by  (35) 

where  k  =  14. 2x  106  [A/to],  and  b  =  7540  [Ns/m\.  The 
initial  simulation  conditions  were  set  to  be  as 

q{t0)  =  0  [C\ ,  Vin  (t0)  =  0  [V]  (36) 

y  (to)  =  3  [//to]  ,  y  (t0)  =  0  (37) 

where  t0  =  0  [sec] .  The  voltage  to  charge  hysteresis 
model  is  implemented  using  the  expressions  in  (2) ,  (3) , 
and  (4),  where  the  following  functions  are  utilized  [8] 

/  =  A  tanh  (£T4)  +  y0Vh  (38) 

g  =  JjLfi-ie-eW)  (39) 

where  <5  =  0.1,  A  =  5000,  £  =  3/5,  yQ  =  0.1,  7  =  3/5, 
(3  =  0.1,  and  tanh  (.)  is  the  standard  hyperbolic  tangent 
function.  The  expressions  selected  for  /  (V/)  and  g  (V/) 
in  (38)  and  (39)  were  selected  to  meet  Properties  1-5. 
The  following  desired  trajectory  was  selected  for  the 
tracking  problem 

yd  =  1.5  x  10'6  sin(27 rft  -  0.5t t)  +  1.5  x  10~6.  (40) 

It  is  clear  from  (40)  that  Assumption  4  is  met.  To 
demonstrate  proof  of  concept  for  this  approach,  three 
frequencies  (10[ilz],  100 [Hz],  and  1000[Hz])  were  se¬ 
lected  in  an  attempt  to  demonstrate  bandwidth  of  the 
proposed  approach.  Table  1  shows  the  applied  de¬ 
sired  frequency  (ADF),  sampling  frequency  (SF),  con¬ 
trol  gains,  and  time  needed  to  catch  up  with  the  desired 
trajectory  (TT).The  PZTA  desired  trajectory  yd(t) 
and  actual  displacement  y  (t),  tracking  error  e  (/),  and 
voltage  V),  ( t )  are  presented  in  Figures  3,  4,  and  5,  re¬ 
spectively,  for  the  selected  frequency  /  =  10  [Hz],  Fig¬ 
ures  6,  7,  and  8,  respectively,  for  the  selected  frequency 
/  =  100  [Hz],  and  Figures  9,  10,  and  11,  respectively, 
for  the  selected  frequency  /  =  1000  [Hz].  It  is  clear 
from  these  results  that  this  control  strategy  meets  the 
desired  objective. 

5  Conclusion 

This  work  develops  a  nonlinear  robust  controller  for  a 
piezoelectric  actuator  where  its  effective  tip  is  driven 
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ADF  [Hz] 

10 

100 

1000 

SF  [Hz] 

io-y 

3xl0~s 

10^Y 

a 

1,000,000 

140,000 

100,000 

kr 

300,000 

100,000 

50,000 

ks 

10-7 

10^5 

100 

£ 

10-7 

HP5 

IP5 

p 

1000 

1000 

100 

TT  [Sec] 

2.64xl0~4 

1.9xl0-4 

0.8xl0-4 

Table  1:  Simulation  parameters  for  different  frequencies 
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Figure  3:  PZTA  desired  trajectory  yd  ( t )  vs.  actual  dis¬ 
placement  y  ( t )  for  a  desired  frequency  of  10 Hz 

to  track  a  desired  trajectory.  The  PZTA  charge  feed¬ 
back  along  with  the  partial  knowledge  of  the  hysteresis 
model  is  utilized  to  design  a  nonlinear  robust  control 
strategy.  A  Lyapunov-based  analysis  is  presented  for 
the  control  design  development.  Simulation  results  are 
presented  to  illustrate  the  proof  of  concept  of  the  pro¬ 
posed  control  scheme. 
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Appendix 

A  Proof  of  Property  6 

Lemma  1:  Based  on  the  definitions  presented  in 
(3)  and  (4),  the  function  d(Vh)  can  be  upper 
bounded  as  shown  in  (5)  where  r/d  £  R+  is  a 

constant  if  lim  d(Vh)  =  lim  d(Vh)  =  0. 

Vh >  —  oo  Vh >oo 

Proof:  Based  on  the  status  of  Vh  (t),  two  cases  arise 
for  the  hysteresis  model  as  in  (3)  and  (4). 


Case  I:  Vh  (f)  is  increasing  and  14  (t)  >  T4o 

When  Vh  (t)  is  increasing  and  Vh  (t)  >  Vflo ,  the  limit  of 
the  expression  defined  in  (3)  can  be  written  as  follows 

lim  ([q0-f(Vh J]  (41) 

Vh - too  \ 

Vh 

+e~SVh  j  ( g  (t)  -  /'  (r))  eSTdr 
vho 

Because  q0  and  /  (I4J  are  constants,  the  limit  in  (41) 
can  be  rewritten  as  follows 

fao-mj]  lim  e-s(-Vh~Vho)  (42) 

Vh — yoo 

Vh 

+  lim  e~SVh  lim  [  ( g  (t)  —  f  (t))  eSTdr. 

Vh  — yoo  Vh— yoo  J\  V  7 

Vho 


5,  it  is  clear  that  the  upper  bound  in  (43)  is  applied. 
Hence,  it  is  clear  that  lim  d  (Vh)  =  0.  From  the 

Vh — y-oo 

results  obtained  in  Cases  I  and  II,  it  is  obvious  that 
| d  (Vh)  |  <  r)d,  where  r/d  €  K+  is  a  constant. 

B  Proof  of  Theorem  1 

Let  V ( z ,  t)  €  R.  denote  the  following  non-negative  func¬ 
tion 

V  =  ^ e 2  +  ^mr2  (46) 

where  z(t)  was  defined  in  (32).  Note  that  (46)  is 
bounded  as  follows 

Xi\\zf  <V(z,t)<X2\\z\\2  (47) 

where  Al5A2  €  R+  are  constants.  After  taking  the 
time  derivative  of  (46) ,  the  following  expression  can  be 
written  as 


It  is  clear  that  the  limits  lim  e  s^Vh  Vh°  1  and 

Vh — yoo 

lim  e~SVh  equal  to  zero.  From  Properties  4  and 

Vh — yoo 

5,  it  is  clear  that 


J  ( 9  (r)  -  /  (r))  eSrdr 
vho 


<  V  i 


(43) 


where  €  R+  is  a  constant.  Hence,  it  is  clear  that 
lim  d(Vfc)  =  0. 


Case  II:  V),  (t)  is  decreasing  and  V),  (t)  <  Vha 


V  =  -ae2-krr2+(^^j  rs+rNd+rN-^p(\\z\\)2  ||2i||2r 

(48) 

where  (12),  and  (34)  were  utilized.  By  utilizing  (33), 
and  the  triangle  inequality,  V  ( t )  can  be  upper  bounded 
as  follows 

V  <  -a\e\2  ~kr\r\2  +  (^^\\r\\s\+\r\\N^9) 

+  kl p(II^II)  \\4  -  7p(IMI)2  INI2  M2- 

Note  that  the  terms  |r|  |s|  and  |r|  |A4|  can  be  upper 
bounded  as  follows 


When  Vh  (t)  is  decreasing  and  Vh  ( t )  <  I4o ,  the  limit  of 
the  expression  defined  in  (4)  can  be  written  as  follows 


MM  <  T-  M2  +  $i  M2  (50) 

oi 

|r|  |JVd|  <  ^\r\2  +  S2\Nd\2 

o2 


lim  ([q0^f(Vho)\e-s^-v^  (44) 

vh — y-oo  \ 

vh 

+eSVh  J  (g  (r)  -  /  (r))  e^dTdr 

Vho 

By  doing  the  same  procedure  done  for  Case  I,  the  limit 
in  (44)  can  be  rewritten  as  follows 

(qo~f(Vh  J]  lim  e-W..-^)  (45) 

Vh — ►- oo 

Vh 

+  lim  eSVh  ^  lim  J  (g  (r)  —  /  (r)^  e~STdr. 

Vh0 


It  is  clear  that  the  limits  lim  e  s(Vh<>  Vh)  and 

Vh — >-oo 

lim  eSVh  equal  to  zero.  From  Properties  4  and 
vh — y  oo 


where  hi, S2  €  R+  are  constants.  By  utilizing  (50), 
V (t)  of  (49)  can  be  upper  bounded  as  follows 

V  <  -A3  N|2+p(N|)  \\z\\  |r|  (l  -  \p(\\z\\)  |M||  \r\)+e, 

(51) 

where  A3=  min  ja,  (jcr  —  }  ,  provided  kr  is 

selected  such  that  kr  >  The  variable 

e0  €  R+  is  defined  as  follows 

£o“  (^)|s|2  +  ^l^|2.  (52) 

The  expression  in  (52)  can  be  upper  bounded  by  the 
following 

£o  <  ^  ^  V2  +  S2r/2Nd  (53) 


where  the  following  bounds  are  utilized 


where  (29)  can  be  written  as  follows 


M  <  Vs  \N4  <  vNd 


(54)  N(y,y,e,r,yd)  =  myd  +  FL  +  (l  -  a2m)  e+amr  (60) 


where  r/s,riNd  £  R+  are  constants,  and  (24)  and  Prop¬ 
erty  6  were  utilized.  From  (51),  two  cases  arise. 

Case  I:  (e  <  p(||z||)  ||z||  |r|) 

The  expression  in  (51)  can  be  written  as  follows 

<  -A3  ll-ll2  +  (55) 

where  the  term  p(||z||)  ||z||  \r\  (l  -  jp(p||)  ||z||  |r|)  in 
(51)  will  be  always  negative. 

Case  II:  (e  >  p (||^||)  INI  M) 

The  expression  in  (51)  can  be  written  as  follows 

R<  -A3|M|2  +  £i  (56) 

where  e4  £  R+  is  defined  as  =  e  +  e0. 

In  both  cases,  (55)  and  (56)  can  be  written  as  follows 

VK-^-V  +  e!.  (57) 

A2 

From  (57)  it  is  clear  that  ||z||  is  upper  bounded  as  fol¬ 
lows 


where  (12)  was  utilized.  To  further  facilitate  the  sub¬ 
sequent  analysis,  the  following  terms  N(y,yd,0,0,yd), 
N(y,  y,  0, 0,  yd),  and  N(y,  y,  e,  0,  yd)  are  added  and  sub¬ 
tracted  to  the  right-hand  side  of  (28)  as  follows 

N  =  [N(y,yd,0,0,yd)  -  Nd(yd,yd,0,0,yd)]  (61) 

+  [N(y,  y,  0, 0,  yd)  ~  N{y,  yd,  0, 0,  yd)\ 

+  [N(y,  y,  e,  0,  yd)  -  N(y,  y,  0, 0,  yd)\ 

+  [■ N(y ,  V,  e,  r,  yd)  -  N(y,  y,  e,  0,  yd)\ . 


After  applying  the  Mean  Value  Theorem  to  each  brack¬ 
eted  term  of  (61),  the  following  expression  can  be  ob¬ 
tained 


N  = 


dN(a1,yd,0,0,yd) 


da  1 

dN(y,a2,0,0,yd) 


+ 


da-2 

dN(y,y,a3,0,yd) 


da3 

dN(y,y,ei,a4,yd) 


da4 


(V  ~  Vd)  (62) 

(y  -  Vd) 

’ 2 

(e-0) 

} 3 

(r-0) 


where  v4  £  {yd,y),  v2  £  (yd,y),  v3  £  (e,  0),  and 
v4  £  (r,  0).  The  right-hand  side  of  (62)  can  be  upper 
bounded  as  follows 


INI  <  VPo  exp  (-Pit)  +  p2  (1  -  exp  (58) 

where  (3a  =  yf  Wz  (*o)||  ,/?i  =  yf,  and  (32  = 

From  (58)  it  is  clear  that  2  (t)  is  practically  regulated 
to  0  as  t  — >  00,  hence,  from  (32)  e(t)  is  practically 
regulated  to  0  as  t  — >  00,  thus  meeting  the  proposed 
control  objective.  To  meet  our  secondary  objective,  it  is 
clear  from  (47),  (55),  and  (56)  it  is  clear  that  V(z,t)  £ 
Coo;  hence  z(t) ,  r  (t)  ,e(t)  £  C^.  From  (12),  it  is  clear 
that  e  (t)  £  Coo-  Thus,  from  (13)  and  Assumption  4,  it 
is  clear  that  y{t),y  (t)  £  Coo-  Since  z[t)  ,r  ( t )  €  Coo, 
then  from  (31),  it  is  clear  that  qd  (t)  £  Coo-  From  the 
bounding  statement  in  (19)  and  (24),  it  is  clear  that 
q(t)  £  Coo,  and  then  Vm{t)  £  Coo-  Thus,  from  (21), 
it  is  clear  that  Vh(t)  £  Coo-  From  (6)  and  (8),  it  is 
clear  that  Vc(t)  £  Coo,  thus,  it  is  clear  that  V  ( t )  £ 
Coo-  From  (1),  and  (7),  it  is  clear  that  y  (t)  £  Coo- 
Standard  signal  chasing  arguments  can  be  used  to  prove 
all  remaining  closed-loop  signals  are  bounded. 

C  Upper  Bound  Development 

Recall  the  definitions  in  (29)  and  (30) 

N  =  N(y,  y,  yd)  =  myd  +  FL+e  +  ame  (59) 
Nd  =  Nd(t)  =  myd  +  F  Ld 


dN(ai,yd,  0, 0,  yd) 

da± 

O’  l=Vl 

dN(y,  a2, 0, 0,  yd) 

da  2 

O2—V2 

dN(y,  y,  a3,0,yd) 

da3 

CT3=V3 

dN(y,  y,  e,  a4,yd) 

da4 

(T4—V4 

(63) 


The  partial  derivatives  in  (63)  can  be  calculated  from 
(60)  as 


dN(a1,yd,0,0,yd ) 
da  1 

=  fL  1 

(64) 

dN(y,  o-2,0,0,  yd) 
da  2 

=  /i2 

(65) 

dN(y,  y,a3, 0,  yd) 
da3 

=  1  —  a2m 

(66) 

dN(y,y,e,a4,yd) 

da4 

=  am 

(67) 

where  (10)  was  substituted  in  (60),  and  the  functions 
/li  (a  1,  yd)  €  R  and  /l2  (y,  a2)  £  R  are  defined  as 
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follows 


f L1 

a  dFL  (a  1 ,  yd) 
da± 

,  Tem 

Cc 

(68) 

f L'2 

a  dFL(y,a2) 

(69) 

da2 

By  noting  that 

vi  =  y-ci{y-yd)  v2  =  y-  c2  {y  -  yd) 

i>3  =  e  —  C3  (e  —  0)  v4  =  v  —  C4  (r  —  0) 

where  Cj  €  (0,1)  V*  =  1,2,. ..,4.  From  Assumptions 
1,  3,  and  4,  it  is  clear  that  an  upper  bound  for  the 
right-hand  side  of  (64)-(67)  can  be  written  as  follows 


dN(a1,yd,0,0,yd) 


da  1 

dN{y,  a2, 0, 0,  yd) 


da  2 

dN(y,y,a3,0,yd) 


da3 

dN(y,y,e,a4,yd) 


da4 


<  Pi  (lehlrl)  (70) 

<  P2  (lel  >  M) 


P3 


<  Pa 


where  pt  Mi  =  1, 2, 4,  are  positive  constants.  After 
substituting  (70)  into  (63),  N(-)  can  be  expressed  as 

N  <  pA\\z\\)\e\+p2(\\z\\)\e\  (71) 

+P3  lel  +  Pa  M  • 


The  expressions  in  (12),  and  (32)  can  now  be  used  to 
upper  bound  the  right-hand  side  of  (71)  as  in  (33). 
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